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I. INTRODUCTION 

Absorbing phase transition (APT) are a particular 
class of non-equilibrium phase transitions occurring in 
physical, biological, as well as chemical systems (see for 
instance 0). Transitions to absorbing states are of par- 
ticular interest since they have no equilibrium counter- 
parts and may occur even in one-dimensional systems. 
A characteristic feature of absorbing phase transitions is 
the competition between the proliferation and annihila- 
tion of a certain entity A, e.g., particles, energy units, 
viruses, molecules in catalytic reactions, etc. It is essen- 
tial that no spontaneous creation of such quantities takes 
place. At a critical value of the proliferation-annihilation 
rate the density p{A) vanishes and the system is trapped 
forever in the absorbing state p{A) = 0. 

Directed percolation is recognized as a paradigmatic 
example of absorbing phase transitions. This is reflected 
by the universality hypothesis of Janssen and Grass- 
berger that models which exhibit a continuous phase 
transition to a single absorbing state generally belong 
to the universality class of directed percolation 0, U 
Different universality classes occur for instance in the 
presence of additional symmetries. In particular, par- 
ticle conservation may lead to the different universality 
class of absorbing phase transitions with a conserved field 
as pointed out in For instance the conserved lattice 
gas (CLG) 0, the conserved threshold transfer process 
(CTTP) (J|, the weU known Manna sandpilc model 0, as 
well as a reaction-diffusion model belong to this uni- 
versality class 0- Note that this universality class is of 
particular interest since the corresponding systems con- 
nect the critical behavior of the absorbing phase transi- 
tion with the critical steady state of self-organized critical 
(SOC) systems Actually, SOC sandpile models can 
be considered as driven-dissipative versions of (closed) 
systems exhibiting absorbing phase transitions |^ . 

In this paper we consider the universal finite-size scal- 
ing as well as the universal dynamical scaling behavior 
of the CLG model, the CTTP, and the Manna model. 



First we introduce a method that allows to study finite- 
size effects in the steady state. In contrast to previ- 
ous attempts to measure finite-size effects our method 
is well defined. Furthermore it can be applied imme- 
diately to other classes of absorbing phase transitions. 
Second we consider the activity spreading of a single ac- 
tive seed. The corresponding spreading exponents are 
naturally connected to the avalanche exponents of SOC 
systems [To| . In particular we discuss certain hyperscal- 
ing laws relating the spreading exponents to the steady 
state exponents of absorbing phase transitions. This al- 
lows us to express the SOC avalanche exponents in terms 
of the exponents of the corresponding absorbing phase 
transition (e.g. the exponents of the order parameter, 
the correlation length exponent, etc.). Thus the criti- 
cal state of SOC systems is closely related to the critical 
properties of an ordinary second order phase transition. 



II. MODELS 

The first considered model is the conserved lattice gas 
(CLG) 1 which is a stochastic variant of a model in- 
troduced by Jensen ^j. In the CLG model lattice sites 
may be empty or occupied by one particle. In order to 
mimic a repulsive interaction a given particle is consid- 
ered as active if at least one of its neighboring sites on the 
lattice is occupied by another particle. If all neighboring 
sites are empty the particle remains inactive. Active par- 
ticles are moved in the next update step to one of their 
empty nearest neighbor sites, selected at random. 

The second model is the so-called conserved threshold 
transfer process (CTTP) 0, a modification of the thresh- 
old transfer process introduced in fl^. Here, lattice sites 
may be empty, occupied by one particle, or occupied by 
two particles. Empty and singly occupied sites are con- 
sidered as inactive whereas double occupied lattice sites 
are considered as active. In the latter case one tries to 
transfer both particles of a given active site to randomly 
chosen empty or singly occupied nearest neighbor sites. 
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It was shown recently that the order parameter and its 
fluctuations obey in the steady state the scaling forms 
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FIG. 1: Sketch of the dynamics of the three considered mod- 
els. Filled circles mark active particles whereas non-active 
particles are marked by open circles. The arrows denote how 
the active particles are (probably) moved in the next update 
step. In the case of the one-dimensional CLG model the par- 
ticle transfer is deterministic. For the one-dimensional CTTP 
stochastic (left) as well as deterministic (right) particle move- 
ments may occur. Only the one-dimensional Manna model is 
characterized by a full stochastic dynamics in D = 1. 



The third model is a modified version of the Manna 
sandpile model |^ the so-called fixed-energy Manna 
model 3- III contrast to the CTTP the Manna model 
allows unlimited particle occupation of lattice sites. We 
use in our investigations the original Manna relaxation 
rules, i.e., lattice sites which are occupied by at least two 
particles are considered as active and all particles are 
moved to the neighboring sites selected at random. 

The three models are sketched in Fig.Q] We use in all 
cases periodic boundary conditions, i.e., closed systems 
are considered and the number of particles is conserved. 
In our simulations (see (I^ for details) we start from 
a random distribution of particles. All models reach after 
a transient regime a steady state which is characterized 
by the average density of active sites p^. The density p^ 
is the order parameter and the particle density p is the 
control parameter of the absorbing phase transition, i.e., 
the order parameter vanishes at the critical density pc 
according to p.^ cx 5p^ , with the reduced control param- 
eter 5p = pI Pc — ^- Additionally to the order parameter 
we consider its fluctuations Ap^. Approaching the tran- 
sition point from above {5p > 0) the fluctuations diverge 
according to Ap^ oc 5p~'' (see Below the criti- 

cal density (in the absorbing state) the order parameter 
as well as its fluctuations are zero in the steady state. 

Similar to equilibrium phase transitions it is possible 
in the case of absorbing phase transitions to apply an 
external field h which is conjugated to the order parame- 
ter, i.e., the field causes a spontaneous creation of active 
particles (see for instance jlj). A realization of the exter- 
nal field for absorbing phase transitions with a conserved 
field was developed in where the external field trig- 
gers movements of inactive particles which may be ac- 
tivated in this way. At the critical density p^ the order 
parameter and its fluctuations scale as p^ oc h^/'^ and 



The universal scaling functions R{x,y) and D{x,y) are 
the same for all systems belonging to a given univer- 
sality class whereas all non-universal system-dependent 
features (e.g. the lattice structure, the update scheme, 
etc.) are contained in the so-called non-universal metric 
factors Op, a^, and |15|. The universal scaling func- 
tions are normed by the conditions i?(l,0) = i?(0, 1) = 
£'(0, 1) = 1 and the non- universal metric factors can be 
determined from the amplitudes of 



Pa(<5p,/i = 0) 
Pa(<5p = 0,/l) 

Oa Ap^{6p = 0, /i) 



(3) 
(4) 
(5) 



These equations are obtained by choosing in the scaling 
forms [Eqs. I2J)] a^SpX — 1 and a^hX"' = 1, respec- 
tively. 



A recently performed analysis of the universal scaling 
functions as well as of the critical exponents for D > 2 Q 
confirms the conjecture of Q that the CLG model, the 



t 



t + 2. 



o o 



o o 



o o 



o o 



o o 



o o 



o o 



o o 



• o 



o o 



o o 



ooo.o.ooooo, oo 



o o 



o o 



o o 



o o 



• o 



FIG. 2: Sketch of the dynamics of the one-dimensional CTTP. 
At update step t an active site is triggered by the external 
field in a cluster of inactive sites. The arrows denote how 
the active particles (full circles) are moved in the next update 
step. Due to the dynamic rules of the CTTP only two dif- 
ferent relaxation processes (selected randomly) occur: both 
particles are moved to the same empty site or both particles 
are moved to two adjacent sites. Thus the dynamics of the 
one- dimensional CTTP is characterized by a trivial random 
walk of the perturbation. This random walk proceeds until it 
reaches the boundary of a cluster. 
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CTTP, and the Manna model belong to the same univer- 
sality class. The situation is more complicated in one- 
dimensional systems where a splitting of the universality 
class occur. The reason for the non-universal behavior 
is that the dimensional reduction changes the stochas- 
tic character of the dynamics (see Fig.H). For instance 
the CLG model is characterized by deterministic toppling 
rules in D = 1 and exhibits a trivial phase transition 
with (3 = 1 and = 1/2 (see also ref. 22 in 0). Due to 
the trivial behavior of the one-dimensional CLG model 
we consider in this work the two- and three-dimensional 
CLG model only. 

In the case of the CTTP we observe that roughly 40% 
of the relaxation events are deterministic. Furthermore 
a perturbation that is triggered by the external field per- 
forms a simple random walk (see Fig.[2Jl. This pathologic 
behavior is completely different from the behavior of the 
one-dimensional Manna model that is characterized by a 
pure stochastic relaxation of active particles to the next 
neighbors. More than the other models the Manna model 
is therefore the paradigm of the universality class of ab- 
sorbing phase transitions with a conserved field. In the 
following we will call that class the Manna universality 
class since universality classes are often labeled by the 
simplest model belonging to them. 

The universality splitting of the one-dimensional sys- 
tems is in full agreement with the universality hypoth- 
esis of sandpile models 0|. According to this conjec- 
ture the universality classes of sandpile models are de- 
termined by the way the particles are distributed to the 
next neighbors (deterministic, stochastic, directed, undi- 
rected, etc.). Obviously the Manna universality class is 
characterized by a stochastic and undirected distribution 
of particles. 



III. STEADY-STATE FINITE-SIZE SCALING 

Similar to equilibrium critical phenomena we assume 
that the system size L enters the scaling forms [Eqs. H1I2|| ] 
as an additional scaling field, i.e. 

p..{Sp,h,L) ^ X-^'R^Ua.SpX^a.hX-' ^a^LX-"^) , 
a^Ap_{dp,h, L) D^i,,{apdpX,a,,hX'' ,aLLX^''^){6) 

where the exponent describes the divergence of the 
spatial correlation length, i.e., oc 6p~'^-^. Note that 
the universal scaling functions depend now on the partic- 
ular choice of the boundary conditions, the system shape 
etc. ^15]. Throughout this work we use in all dimen- 
sions hyper cubic lattices with periodic boundary con- 
ditions (pbc). However, the universal scaling functions 
[Eqs. H1I2|| ] are recovered in the thermodynamic limit, e.g. 



defined as (see for instance [T3) 



(7) 



Q = 1 - 



(8) 



For non-vanishing order-parameter the cumulant tends 
to Q = 2/3 in the thermodynamic limit. In the case of 
a zero order parameter the cumulant vanishes if the or- 
der parameter is characterized by a Gaussian distribution 
symmetrically distributed around zero. The latter case 
is observed in equilibrium systems, e.g. the Ising model 
for T > T„. In the case of absorbing phase transitions 
the order parameter is non-negative per definition. Thus 
the order parameter is characterized by a non-trivial dis- 
tribution close to criticality and the above scenario does 
not apply. 

Nevertheless one expects that the cumulant obeys the 
scaling form 

Q(Jp,/i,L) ^Qp„,(a>A,a,/iA",a^LA-'^^). (9) 

Notice that no metric factor Oq is needed since the cu- 
mulant is already dimensionless. Choosing a^L A^''^ = 1 
we get for zero-field 



Q(0,0,L) = Q{Sp,0,L) 

Sp = 

~ Qp..(a,<5p(a^L)-''^,0,l) 
= Qp..(0,0,l) 



Sp=0 



(10) 



Additionally to the order parameter and its fluctua- 
tions we consider the fourth-order cumulant Q which is 



which is obviously universal. The universal value 
Qpbc(0, 0,1) corresponds to an intersection point if one 
plots Q as a function of p for various system sizes L. 
Thus it is possible to determine the critical value p^ from 
the common intersection point. This cumulant intersec- 
tion method is very useful and was applied in numerous 
works (see for instance ^3 ^^'^ references therein). 

As usually finite-size effects have to be taken into ac- 
count if the correlation length is of the order of the system 
size. A feature of these finite-size effects is that a given 
system may pass within the simulations from one phase 
to the other. This behavior is caused by critical fluctu- 
ations, i.e., approaching the transition point the order 
parameter vanishes whereas its fluctuations diverge. But 
in contrast to common second order phase transitions the 
situation is drastically different in the case of absorbing 
phase transitions. Approaching the transition point the 
correlation length increases and as soon as is of the 
order of L the system may pass to an absorbing state and 
is trapped forever. Additionally to the absorbing phase 
((Pa) = for (5p < 0) the steady state order parame- 
ter and its higher moments vanish {{p^)l — 0) even in a 
small vicinity above the critical point. Thus in the case 
of absorbing phase transitions it is impossible to con- 
sider finite-size effects of steady state quantities around 
the critical point. 

In order to bypass this problem it was suggested (see 
for instance [Tg) to consider metastable (ms) or quasis- 
teady state values of the order parameter and its higher 
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FIG. 3: The decay of the order parameter close to criticahty 
for the two-dimensional CTTP. After a so-called metastable 
regime the system passes to the absorbing phase. The inset 
displays that no clear saturation of the order parameter could 
be observed. 
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FIG. 4: The density of active sites pa as a function of the 
external field for various system sizes L. The inset displays 
the corresponding order parameter fluctuations. The dashed- 
lines correspond to the asymptotic behavior of the infinite 
systems [Eqs.|g]|Sll]. 



moments (pf^)i,m3, respectively. This is shown in Fig.|3|for 
the order parameter of the two-dimensional CTTP close 
to the transition point. After a short transient regime the 
system reaches a metastable state where it can spend a 
certain time until it finally enters an absorbing state. In 
the metastable phase the order parameter is expected 
to fluctuate around a well defined average value (/Oa)i,,„s 
which is used for the finite-size scaling analysis. This 
method was applied in previous works (e.g. |3iSll3 and 
the results sound mostly valid. 

Nevertheless this method can be questioned. First, 
there is no well defined average value of the order pa- 
rameter in the metastable regime. This can be seen in 
the inset of Fig.|31 where we scrutinize the data. No clear 
saturation of the order parameter can be observed. Sec- 
ond the method is quite inefficient. In the case of the data 
presented in Fig.|31we used 5 10^ different initial configu- 
rations for L = 32 to get a sufficiently averaged estimate 
of the order parameter. Although roug hly 10" lattice 
updates k, 2000) were performed no clear satura- 

tion could be observed. Thus reliable data for larger 
lattice sizes (L = 64, 128, 256, ...), which are required for 
an appropriate finite-size scaling analysis, can not be ob- 
tained within moderate computer times. Third and fi- 
nally no rigorous prove exists that the metastable order 
parameter moments (pj;}i,m5 scale in the same way as the 
corresponding steady state order parameter moments. 

In contrast to the consideration of metastable phases 
we choose a different method in order to study finite-size 
effects in the steady state. In our simulations we mea- 
sure the order parameter at the critical density {Sp = 0) 
as a function of the conjugated field h for various sys- 
tem sizes. Due to the external field the system can not 
be trapped forever in the absorbing phase. Therefore 



steady state quantities are available for all densities. In 
Fig.0]we present the order parameter and its fiuctuations 
for the two-dimensional CTTP model. The finite-size ef- 
fects, i.e., the deviations from the behavior of the "infi- 
nite" system (L ^ ) can be clearly seen. Note that the 
data for L = 32 are averaged over 3 10^ lattice updates 
and we obtain smooth curves of the order parameter, its 
fiuctuations as well as of the cumulant. Thus the numer- 
ical effort of this method is significantly smaller than the 
above discussed analysis of the metastable regime. 

According to the above scaling laws [Eqs. H6l9|l ] the 
finite-size scaling forms are given by 

Q{0,h,L) ~ Q^UO,a^h{a^Lr/''^,l). (11) 

For the sake of convenience we norm the universal scaling 
function Q^^c by the condition 

Qp.c(0,l,l) = 0. (12) 

Since the metric factor is known from previous simu- 
lations (via Eq. (Q) the above condition can be used 
to determine the metric factor a^. Taking into account 
that the correlation length scales at criticality as 

a^e^ ~ {a.h)-''^/^ (13) 

we find that Eq. (|12|l implies that the universal func- 
tion Qpb^ is positive for a^L > a±^± and negative for 
a^L < . Note that in the case of equilibrium phase 

transitions Eq. H12|) is useless since the cumulant is usu- 
ally positive. 
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FIG. 5: The universal finite-size scaling analysis of the or- 
der parameter pa, of the fiuctuations Apa as well as of the 
cumulant Q for the one dimensional CTTP and the one- 
dimensional Manna model. The long-dashed lines correspond 
to the power law behaviors of the infinite system [Eqs.(0]|SJ] 
and to the cumulant limit 2/3, respectively. The data are ob- 
tained from simulations of system sizes L G {2048, 4096, 8192} 
where up to 10^'' lattice update steps are performed. 

In Figs . 15171 we present the universal finite-size scaling 
analysis of the CLG model, the Manna model and the 
CTTP for D = 1,2,3. Since the exponents^ and a 
were already determined in previous works 0, we 
just vary the correlation length exponent in order to 
produce data collapses. The value of the non-universal 
metric factor Ql is determined via Eq. 112|l . We observe 
good data collapses for i/j^ — 0.799 ±0.014 for D — 2 and 
= 0.593 ±0.013 for D = 3, respectively. The data col- 
lapses are quite sensitive for variations of the exponents. 
Thus the quality of the corresponding data collapses are 
used to estimate the error-bars. The values of the expo- 
nents as well as of the non-universal metric factors are 
listed in TableHland TableHil 

In the case of the one dimensional Manna model and 
the one dimensional CTTP we observe the expected split- 
ting of the universality class. The correlations length ex- 
ponent Vj^, the field exponent a, as well as the scaling 
functions differ clearly (see Tableland Fig.[Sl). Further- 
more the value of the Manna model i'^ = 1.347 ± 0.091 
differs clearly from t/^ = 1.80 ± 0.01 obtained in a 
previous work including a finite-size scaling analysis of 
metastable states [l9| . 

Despite this splitting of universality we observe that 
the fourth-order cumulant tends for all models in all 
dimensions to infinity if one approaches the transition 
point, i.e., 

Qpbc(0,x, 1) ^ -oo for a; 0. (14) 

This behavior is caused by the vanishing steady state 
fluctuations [see Eq. ©]. Thus we assume that the di- 



FIG. 6; The universal finite-size scaling analysis of the 
order parameter p^, of the fluctuations Ap^ as well as 
of the cumulant Q for the two dimensional models. The 
long-dashed lines correspond to the power law behaviors of 
the infinite system [Eqs.lQ] |SJ] a^nd to the cumulant limit 
2/3, respectively. The data are obtained from simulations 
of system sizes L £ {64, 128, 256} where up to 5 10^ lattice 
update steps are performed. 




FIG. 7: The universal finite-size scaling analysis of the or- 
der parameter p^, of the fiuctuations Ap^ as well as of the 
cumulant Q for the three dimensional models. The long- 
dashed lines correspond to the power law behaviors of the 
infinite system [Eqs.l|l||^] and to the cumulant limit 2/3, re- 
spectively. The data are obtained from simulations of system 
sizes L G {16,32,64} where up to 510^ lattice update steps 
are performed. 



vergent fourth-order cumulant is a characteristic feature 
of all absorbing phase transitions, independently of the 
considered lattice structure as well as of the particular 
considered universality class. Preliminary simulations for 
directed percolation support this conjecture and will be 
published elsewhere. 
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IV. DYNAMICAL SCALING BEHAVIOR 
A. Homogeneous particle source 

In the following we investigate the dynamical scaling 
behavior in the vicinity of the absorbing phase transition. 
First we consider how the order parameter decays, 
starting the simulations from a random distribution of 
particles (so-called homogeneous particle source). Above 
the transition point the density of active sites decreases 
in time and tends to the steady state value (despite of 
finite-size effects as discussed above). Below the transi- 
tion point the density of active sites decreases exponen- 
tially to zero. At the critical point the order parameter 
decays algebraically according to 

P. - (a*t)"" (15) 

where at denotes a corresponding non-universal metric 
factor. A finite system size limits this power-law behav- 
ior and one expects that the order parameter obeys at 
criticality the scaling ansatz 

p^{L,t) ~ A-"--" R'^^MttX-^Ka^LX-^^) (16) 

where we have to distinguish the universal functions R' 
and R. For the sake of simplicity we choose i?'(l, oo) = 1. 
Setting aLLX~'^^ = 1 one gets the finite-size scaling form 

p„(L,<) ^ L-"^i^;,^(a,^(a,L)-^l) (17) 

where z = v^jv^ denotes the dynamical exponent as 
usual. Finite-size effects have to be taken into account 
for 0{t)^t 

Fss with 

tpss = a-^ {a^LY . (18) 
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For t ^ t-psa the scaling function obeys the power-law 
_Rpj^^(x, 1) a;~", whereas R'^^^^{x,l) decays exponen- 
tially for a; ^ 1, i.e., t » tpss- 

According to the above scaling form [Eq. (|17|) ] we plot 
in the Figs. lSllOl the rescaled order parameter as a func- 
tion of the rescaled time. We observe good data col- 
lapses for a = 0.419 ± 0.015, z = 1.533 ± 0.024 for 
D = 2 and a = 0.745 ± 0.017, z = 1.823 ± 0.023 for 
-D = 3, respectively. These values are in agreement with 




a; (a^LY^ 

FIG. 9: The dynamical scaling analysis for the two- 
dimensional models. The long-dashed lines correspond to 
the power law behaviors of the infinite systems [Eqs. 1151 1211 
1^ 1. System sizes from L = 64 up to L = 512 are considered 
and the data are averaged over at least 2 10^ different initial 
natural configurations (see text). 
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FIG. 8: The dynamical scaling analysis for the one- 
dimensional CTTP and the Manna model. The long-dashed 
lines correspond to the power law behaviors of the infinite 
systems [Eqs. (|T3 IHl HU] • System sizes from L = 512 up to 
L — 8192 are considered and the data are averaged over at 
least 10"" different initial natural configurations (see text). 



FIG. 10: The dynamical scaling analysis for the three- 
dimensional models. The long-dashed lines correspond to the 
power law behaviors of the infinite systems [Eqs. ()15ll21ll^^ ]. 
System sizes from L = 16 up to L = 128 are considered and 
the data are averaged over at least 5 10® different initial nat- 
ural configurations (see text). 
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those of previous simulations "4' . In the case of the one- 
dimensional models we observe that the CTTP and the 
Manna model are both characterized by a = 0.141±0.024 
and z = 1.393 ± 0.037 but the corresponding scaling 
curves of both models differ slightly. It is possible that 
this indicates a universality splitting similar to the steady 
state scaling behavior. 

B. Localized particle source 

Additionally to a homogeneously distributed source of 
active sites one usually considers the activity spreading 
generated from a single active seed [l^. In this case it 
is customary to examine the survival probability Pjfy 
that the system is still active after t update steps. Fur- 
thermore one investigates how the number of active sites 
N^{t) increases in time. At criticality the survival prob- 
ability as well as the average number of active sites are 
expected to scale as 

a^N^{L,t) ^ X""^^ N^^^ia,tX-''i^,a,LX-''^). (20) 

The universal functions are normed by P{1, oo) — 1 and 
A'^(l, oo) = 1 and we find in the thermodynamic limit 

apP. - {a,t)-\ (21) 
a^N^ ~ {a,tf. (22) 

The finite-size scaling forms are obtained by setting 
aLLX~^^ = 1, yielding 

a^P.^{L,t) ^ L-^^P,,,(a.^(a^i)-^l), (23) 
a^N^{L,t) ^ L<''N^Uatt{a^L)-\l). (24) 

Again the scaling functions P^^^ and N^^^ decay expo- 
nentially for t ^ ipss whereas they exhibit an algebraic 
behavior for t <^ tpss • 

Since the absorbing state is non-trivial one has to in- 
vestigate the spreading activity at the so-called natural 
density (see for instance J^])- For each considered model 
an absorbing state at is prepared and a particle is 
moved to a randomly selected site in order to create one 
active seed. We obtain for all dimensions convincing data 
collapses which are shown in the insets of the Figs. l8llOI 
The values of the exponents agrees with those of previous 
works Q and are listed together with the non-universal 
metric factors in Table|l] and TablelTll In summary ac- 
tivity spreading from a localized seed is characterized for 
all three models by the same universal scaling functions 
P and TV. 

The activity spreading of APT with a conserved field 
is closely connected to avalanche processes in SOC sys- 
tems 0. In particular the Manna model is a paradig- 
matic example of a class of SOC systems, the so-called 
sandpile models. The SOC version and the APT ver- 
sion of the Manna model are characterized by the same 



TABLE I: The critical exponents of the considered models 
below the upper critical dimension = 4. The data of the 
exponents /? and a are obtained from fT4| . The finite-size 
scaling analysis yields the values of and 7', whereas the 
exponents a, 5, 9, and z are obtained from activity spreading 
(see text). The values of (3' and i'± are determined via scal- 
ing laws. In particular the values of z/y are in good agreement 
with those of direct measurements of the order parameter per- 
sistence distribution [27ll . In the case of the one-dimensional 
models we observe a splitting of the universality class. 





D = 1 


D = 2 


D = 3 


a 


0.141 ±0.024 


0.419 ±0.015 


0.745 ± 0.017 


5 


0.170 ±0.025 


0.510 ±0.020 


0.765 ± 0.025 


9 


0.350 ± 0.030 


0.310 ±0.030 


0.140 ± 0.030 


z 


1.393 ±0.037 


1.533 ± 0.024 


1.823 ±0.023 


13 


0.382 ±0.019 


0.639 ± 0.009 


0.840 ±0.012 


P' 


0.319 ± 0.052Ma„„a 


0.624 ± 0.029 


0.827 ± 0.034 


a 


2.710 ± 0.040Ma„„a 


2.229 ± 0.032 


2.069 ± 0.043 




1.770 ± 0.058cTTP 








1.347 ± 0.09lMa„„a 


0.799 ± 0.014 


0.593 ±0.013 




1.760 ± O.OeOcTTP 








1.876 ±0.135Ma„na 


1.225 ±0.029 


1.081 ±0.027 




2.452 ±0.106cTTP 






i 


0.550 ± 0.040Ma„„a 


0.367 ±0.019 


0.152 ±0.017 




0.670 ± 0.040cTTP 







(microscopic) dynamic rules but the boundary conditions 
differ. Closed boundary conditions lead to a globally con- 
served particle density in the case of absorbing phase 
transitions. Whereas sandpile models are per definition 
driven dissipative systems where particles (sand-grains) 
are injected into the system and dissipated through open 
boundaries. The self-organization of SOC systems cor- 
responds to the fact that they approach, without any 
external fine tuning, the critical state (p(t) — > Pc) in the 
infinitesimally slow driving limit (so-called separation of 
time scales, see [IJ). In the critical state the external 
driving triggers (scale invariant) avalanche- like relaxation 
events. These avalanche processes are described by cer- 
tain critical exponents which can be derived from the 
spreading exponents S, 6, and z 0,^^. In particular the 
avalanches are characterized by several quantities (see for 
instance [l^ |2^'). e.g. the size s (number of elementary 
relaxation events) , the area a (number of distinct toppled 
sites), the time t (number of parallel updates until the 
configuration is stable), as well as the radius exponent r 
(radius of gyration). In the critical steady state the cor- 
responding probability distributions decay algebraically 

P^ oc x-^- , (25) 

characterized by the avalanche exponents with x G 
{s,a, t,r}. Assuming that the size, area, etc. scale as a 
power of each other, 

X oc x''^'^"' (26) 
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FIG. 11: The avalanche exponents of the Manna model in var- 
ious dimensions. In order to avoid overlaps the exponents are 
slightly shifted. The avalanche exponents of the SOC version 
of the Manna model (left) are obtained from [2§] for D — 1 
and from D = 2, 3. Using the Eqs. H29I33^ we ob- 

tained the avalanche exponents (middle) from the spreading 
exponents S, 6, and z. Using certain hyperscaling relations it 
is possible to express the avalanche exponents (right) in terms 
of the exponents of the contirmous absorbing phase transition, 
see Eqs. I54I57II . 



one obtains the scaling relations 

T^, - 1 



7x 



1 



(27) 



The exponent equals the dynamical exponent z, the 
exponent 7^^ corresponds to the fractal dimension of the 
avalanches and the exponent 7^,^ indicates whether mul- 
tiple toppling events are relevant {■y^a > 1) or irrelevant 
= 1). 

These avalanche exponents are connected to the 
spreading exponents S, 9, and z (see for instance |lCl(). 
First the survival probability P^it) is simply given by the 
integrated avalanche duration 



PM = p^t') 



(28) 



yielding 

r, = 1 + 5. (29) 

Since 7^^ = z the radius exponent is given by 

= I + z5. (30) 

Taking into account that the avalanches of the Manna 
model are compact (7 = D) below the upper critical 
dimension A 4 UlEi we find 



1 



z5 



(31) 



Finally the number of topplings St for an avalanche that 
is active at time t equals the integrated numbers of active 



(32) 



t'=0 



leading to |lf 



1 + 61 + 5 



(33) 



Thus the avalanche exponents of the Manna model are 
naturally related to the spreading exponents of the ab- 
sorbing phase transition. In Fig. llll we compare our re- 
sults with those of SOC simulations of the Manna model. 
The data show that the above scaling relations [Eqs. H29t - 
ESJ] are fulfilled. 



V. DISCUSSION 

In this section we check several scaling relations. 
Due to the pathological behavior of the one-dimensional 
CTTP we use the corresponding values of the Manna 
model for our consideration. At the beginning we check 
the scaling relation 



7 



2/3 



(34) 



which can be easily derived from the scaling form of the 
order parameter histogram (see for instance |0). The 
corresponding data are plotted in Fig.^1 As can be seen 
the above scaling relation is fulfilled within the error- 
bars. 
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FIG. 12: Test of the scaling relations 7' = Du±~2i3, a = P/u^ 
(inset), as well as S = P/i^w (inset). 
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Taking into consideration that a weak external field 
may trigger spreading events one finds that the field ex- 
ponent is given by [ij 



(35) 



In Fig. ^1 we check this scaling law. As can be seen it is 
fulfilled within the error-bars. 

Next we consider the hyperscaling relation 



z v., 



(36) 



This scaling relation can be derived if one assumes that 
the steady state scaling forms and the dynamical scaling 
forms can be combined to 

N^{5p, h, L, t, p,^,o) = P^{Sp, h, L, t, p,,o) (37) 

Here the initial density of active sites appears as an 
additional scaling field (see for instance [l| and references 
therein) and the scaling function behaves asymptotically 
as 



^pb,(0,0, CXD, 1,2:) 



X for X 1 
const for a; 1 , 



(38) 



Starting at criticality from a low density of active sites 
(e.g. several seeds) the number of active sites increases as 
iV^ oc until it reaches a maximum and crosses over to 
the expected asymptotic decay oc t~°' . The crossover 
time is determined by 



(39) 
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FIG. 13; Test of the hyperscaling relations g — Du±+V\\ — ^^ii^- 
and 6 — D/ z — — 5 (see inset). 



which corresponds to a merging of the survived (and for- 
mer separated) clusters of activity . The scaling of the 
crossover time explains the choice of the scaling expo- 
nents in Eq. H37|) . Setting atiA~''ii = 1 we find at critical- 
ity 



iV,(0,0,L,i,p.,o)-L^ (a.i)-^/'^" 

i?p,.(0, 0, a^L{a,ty^'', 1, aoPa.o(a*<)^/""*) ■ 



(40) 



The full crossover can be observed if the particular value 
of the initial density of active sites p^ o leads to 1 ^ 
t^o <g; tpss (as it was for instance observed for directed 
percolation [3|). Therefore, the first regime [Eq. (|22|) ] is 
obtained for 1 ^ t ^ the second regime [Eq. itT^ ] 
for t^a <^ tpss and finite-size effects take place in the 
third regime for t > ipss- 

In the case that one starts the simulations with a single 
seed /9,,o = we get t,, oc L^^/(^-^'5). Taking into 

account that z < zD/{D — zS) we find that ipss < ^co, i-e., 
finite size effects take place before the algebraic decay of 
active particles starts [Eq. H15|) ] and the second scaling 
regime does not occur. Furthermore we can use for t <^ 
tpss the approximation 

i?p,.(0, 0, a,L(a,^)-l/^ 1, a„L-'' {a.t)''^'-') 
« i?p,,(0,0,«),l,aoP.,o(M)''/""^) 



T -D I M\D/z 



and Eq. reads 



(41) 



(42) 



Comparing this result with Eq. (|22|) we obtain the hy- 
perscaling relation Eq. H36() as well as Oq = l/a„. In the 
inset of Fig.^l we display the data of the corresponding 
exponents. The hyperscaling relation Eq. H3()|l is fulfilled 
within the error-bars. 

The situation is completely different if one starts the 
simulations with a homogeneous particle source. For 
instance a random distribution of particles leads for 
the two-dimensional CTTP to an initial density p^ g w 
0.1703. In that case the crossover time 



1 f 1 

tco = — Pa.O 

at Va„ 



1.14 



(43) 



is too small and the short time scaling regime {N^ oc 
for 1 ^ < ^ i^„) can not be observed. On the other hand 
the scaling form Eq. 1)37(1 yields for the second regime 

(tco < t < Uss) 

p^{Sp,h,L,t,p,,.o) (44) 
~ (a.t)-^/-^" i?,,.(0, 0, a,L(a,^)'l/^ 1, aoP^.o*""/^"') 
« (a,t)-^/'^iii?p,, (0,0,00,1,00). 

Comparing this result with Eq. H15|l we get the scaling 
relation 



(45) 
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as well as i?pbc(0, 0, cxd, 1, oo) = 1. But as can be seen 
from the inset of Fig.^] this scaling relation is clearly 
violated in _D = 1 and _D = 2. For D — 3 we think 
that the violation of Eq. (|45|l is hidden by the overlap- 
ping error-bars, i.e., the above scaling relation is violated 
below the upper critical dimension, as already observed 
in previous simulations 'i^ . Furthermore the violation of 
the scaling relation Eq. H45|) explains why the well known 
hyperscaling relation of directed percolation 



D 

z 



(46) 



is not fulfilled for absorbing phase transitions with a con- 
served field. 

It is worth to mention that this scaling anomaly is 
not caused by the particular initial configuration (ran- 
dom distribution of particles). We have observed the 
same behavior for a more natural initial configuration, 
where the correlations of the active and non-active sites 
are not trivial. In this case we start the simulations from 
a steady state at the critical density with non-zero field h. 
Switching off the external field we have measured the re- 
laxation of the order parameter from the initial density 
Pa,o = Pa(Pc) h). For p^ Q « 0.1 we observe the same scal- 
ing function and the same exponent a as in the case of a 
random initial configuration. 

A possible explanation of the above scaling anomaly is 
that the asymptotic scaling regime is so far not observed. 
The situation is sketched in Fig.^l In the case of the 
two-dimensional CTTP the crossover takes place at t^o ~ 
1.14 whereas finite-size effects occur at ipss ~ 6000 for 
L ~ 512. Thus the condition <^ t ipss seems to be 
fulfilled and one would expect to observe the asymptotic 



scaling behavior t~^/^\ i . But it is known that crossovers 
span several decades of magnitudes, usually 5, 6, 7 or 
even more decades (see for instance [2^ for a recent work 
on crossover effects for the same universality class). In 
this way it is possible that one observes in simulations a 
smaller exponent {a — 0.419) than the asymptotic value 
(/3/i^ii = 0.522). Future work on larger lattice sizes [L ^ 
512) are needed to clarify whether the scaling anomaly 
can be explained by a simple crossover effect. 

Finally we consider the percolation probability Pperc 
that a path of active sites propagates through the system. 
Obviously the percolation probability is related to the 
survival probability via 



lim P,(t) , 

t — >-oo 



(47) 



The percolation probability is expected to vanish at the 
critical density according to 



(48) 



Assuming that the survival probability obeys the scaling 
form 

apP^{Sp,L,t) (49) 
- A-^^^ii P.U^.Sp^a.tX-'^a^LX-''^) 
we find in the thermodynamic limit 
Pp„o = lim P,{Sp,t, oo) 

~ a;i(a»^'^ii Pp,.(l,oo,oo) (50) 
leading to the scaling relation 

/3' 



S = 



(51) 













N~t \ \\ 







In t 



FIG. 14: Sketch of the crossover behavior of the number of 
active sites A'^^. For t <C tco scales as whereas it is 
expected to decrease for t S> too as cx . It is possible 

that a too low exponent (a < l3/v\\) is observed in simulations 
if one has not reached the asymptotic scaling regime (long- 
dashed line). 



Thus the hyperscaling relations Eq. (|35l I36|l read now 

a = u^D + j^ii - 13' , (52) 



v^D - f3 - f3' 



(53) 



In the case of directed percolation (3 equals /3' due to a 
special symmetry under time reversal (2Q|] . In a more gen- 
eral context one expects that both exponents differ 0], 
for instance in systems with infinitely many absorbing 
states or as in our case, in systems where a conserved 
field couples to the order parameter. The number of in- 
dependent critical exponents is therefore expected to be 
four (e.g. /3, i^y, v^) instead of three independent ex- 
ponents for directed percolation (/3, i/y, Vj^). In order 
to check this scenario we compare in the inset of Fig.ll2l 
the spreading exponent 5 with Surprisingly we ob- 

serve that both values agree within the error-bars for all 
dimensions, suggesting (3 = (3' . It is possible that the un- 
certainty of our estimates hides a tiny difference of both 
exponents. A more accurate test of the scaling relation 
(3 = (3' could be obtained by a direct measurement of the 
percolation probability which remains the topic of future 
research. 
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TABLE II: The non-universal metric factors of the considered 
models. The uncertainty of the metric factors is less than 5%. 





D = 1 


D — 2 


D = 3 


CTTP 








«P 


0.607 


0.341 


0.384 


ah 


0.220 


0.013 


0.093 




187.7 


45.42 


24.51 


at 


1014. 


4.617 


2.173 


at 


1379. 


24.90 


4.239 


ap 


0.107 


0.078 


0.094 


On 


6.062 


2.818 


1.069 


CLG 








ftp 




0.509 


0.434 


ah 




0.062 


0.391 






9.241 


8.881 


at 




2.107 


1.441 


at 




11.22 


3.140 


ap 




0.157 


0.183 


Ojv 




1.249 


0.569 


Manna 








Op 


0.662 


0.211 


0.311 


ah 


7.52 10"^ 


0.007 


0.074 




588.9 


78.56 


32.24 


at 


205.9 


6.011 


2.367 


at 


7.99 10* 


35.53 


4.824 


ap 


0.063 


0.059 


0.089 


aN 


64.02 


3.600 


1.229 



VI. CONCLUSIONS 

We analyzed numerically the critical behavior of three 
different models exhibiting a continuous phase transition 
into an absorbing state. In particular we introduce a 
method which allows to consider finite-size effects in the 
steady state. It is therefore possible to obtain accu- 
rate estimates of the correlation length exponent. Ad- 
ditionally we determine the spreading exponents which 
describe the spreading of activity at the critical point. A 
detailed analysis of certain scaling relations shows that 
usual hyperscaling relations are fulfilled. Only the ac- 
tivity spreading from a homogeneous particle source ex- 
hibits a scaling anomaly. So far this scaling anomaly is 
not understood and remains the topic of future research. 
The number of independent exponents is due to the scal- 



ing anomaly of the exponent a at least four. In the case 
that (3 ^ P' it is even five. 

Since the hyperscaling relations Eqs. (|35l I36|l are ful- 
filled it is possible to connect the SOC avalanche expo- 
nents to the steady state exponents of the corresponding 
absorbing phase transition: 



Tt = 



1+ z +D 

D a- 
2+ 

z Z^ii 

1 + ^^. 

2+ — - 



D Dv^ 
z^ii -I- v^D - a 
i^ii + ly^D - P 



(54) 



(55) 



(56) 



(57) 



In Fig.^] we compare these values with the avalanche 
exponents obtained from SOC simulations of the Manna 
model mmiiil. AU SOC exponents agree within the 
error-bars with the avalanche exponents derived via the 
above scaling laws. Thus it is possible to express the 
avalanche exponents {t^, t„,. . .) of SOC systems in terms 
of the usual critical exponents of a second order phase 
transition (/3, i^y, . . .). In this way, the critical state 
of SOC sandpile models is closely related to the critical 
state of an ordinary second order phase transition. 
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